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Abstract 

Light dark matter (~ 1 — 30 MeV) particles pair produced in electron-positron annihilation 
e“e+ XX inside the supernova core can take away the energy released in the supernova SN1987A 
explosion. Working within the formalism of g-deformed statistics (with the average value of the 
supernovae core temperature (fluctuating) being Tsn = 30 MeV) and using the Raffelt’s criterion 
on the emissivity for any new channel e(e"’'e“ —)> xx) ^ 10^® erg we find that as the 

deformation parameter q changes from 1.0 (undeformed scenario) to l.l(deformed scenario), the 
lower bound on the scale A of the dark matter effective theory varies from 3.3 x 10® TeV to 3.2 x 10^ 
TeV for a dark matter fermion of mass m^ = 30 MeV. Using the optical depth criteria on the free 
streaming of the dark matter fermion, we find the lower bound on A 10® TeV for m^^ = 30 MeV. 
In a scenerio,where the dark matter fermions are pair produced in the outermost sector of the 
supernova core (with radius O.Oiic ^ r < Re, Rc{= 10 km) being the supernova core radius or 
the radius of proto-neutron star), we hnd that the bound on A (~ 3 X 10^ TeV) obtained from 
SN cooling criteria (Raffelt’s criteria) is comparable with the bound obtained from free streaming 
(optical depth criterion) for light fermion dark matter of mass m^- = 10 — 30 MeV. 
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I. INTRODUCTION 


The fact that dark matter plays a very important role in building the universe we live in 
is gradually gaining ground due to concrete experimental evidence collected over a period of 
time. In 1933, Fritz Zwicky [1] found that the normal luminous matter is alone not sufficient 
to explain the velocity dispersion of galaxies in the Coma cluster of galaxies: one requires 
non-luminous matter, dubbed as dark matter (DM). Current data suggests that DM is hve 
times more than the normal luminous matter in the Universe[2]. Now, the dark matter(DM) 
does not interact electromagnetically with the normal luminous matter since it (DM) has no 
electromagnetic charge. Even it does so, it is very weak. So far, scientists have been able to 
infer the existence of dark matter only through its gravitational effect on normal matter. 

But what is dark matter? For a long time it remains a mystery. A wide range of collider 
and astrophysical study suggests that it is a Weakly Interacting Massive ParticIe(WIMP) of 
mass lying in between few MeV to few GeV. Theories suggest that DM candidates are most 
likely to be found in the beyond the Standard Model(SM) physics, such as supersymmetry 
and extra dimensions. Direct detection of DM includes its interaction with nucleons in 
underground detectors, whereas indirect detection through DM annihilation to SM states 
in the Sun (to neutrinos) has been done. Experiments at the Large Hadron Collider(LHC) 
and the upcoming electron-positron linear collider(LC) will give more information about the 
dark matter as the missing energy signature. See HU for a review on dark matter searches. 

An enormous amount of gravitational binding energy 10^^ ergs was released in the su¬ 
pernova SN1987A explosion of which about 99% was carried away by neutrino alone. To 
understand the supernova energy loss mechanism and its relevance in the beyond standard 
model physics has been an area of active research for a long time|5]. In 2006, Fayet et a/. |6] 
studied the impact of light dark matter on the core collapse supernova cooling and found 
that the 1 — 30 MeV mass dark matter fermion can explain SN1987A energy loss rate. They 
also found that if the dark matter particles are of mass < 10 MeV and reproduce the 
observed dark matter relic density, it would lead to the modification of the supernova cool¬ 
ing dynamics, which is unacceptable. Kadota et al. na studied the impact of International 
Linear Collider(ILC) and SN1987A energy loss rate due to the light (MeV) dark matter. 
They found the SN bound to be more stringent than those obtained from ILC by a factor 
of (P(10®) for a DM mass below 100 MeV. 
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Here we would like to investigate the impact of light dark matter fermions on the energy 
released in SN1987A explosion. We work in a dark matter model characterized by an effective 
scale A and use the formalism of g-deformed statistics which takes care of the fluctuation of 
the core temperature of the supernova. While Kadota et al. m considered the dark matter 
fermion and SM photon coupling arising from magnetic dipole moment operator, we have 
generalized our work where a SM photon may couple with dark matter fermion through the 
magnetic dipole moment operator or electric dipole moment operators or both. 

The outline of the work is as follows. We give a brief description of SN1987A cooling 
problem and a small introduction of g-defomed statistics in section II. In section III, we 
discuss the dark matter pair production in electron-positron annihilation and hnd the su¬ 
pernova energy loss rate due to this dark matter fermion pair production. The numerical 
analysis part is presented in section IV. Using the Raffelt’s criterion, we obtain bound on 
the scale A of the effective dark matter theory in deformed(g 7 ^ 1) and undeformed (g = 1) 
scenarios. Using the optical depth criteria(based on free streaming of dark matter fermions), 
we obtain constraints on A. Finally, we conclude in section V and VI. 

II. SUPERNOVA EXPLOSION, ITS COOLING AND g-DEFORMED STATISTICS 
A. Supernova cooling, Raffelt’s criterion 

The supernova SN1987A, a typical example of a core-collapse supernova explosion is the 
hnal fate of massive star of mass M > 8Mq. The energy released in SN1987A explosion 
is enormous: it is the gravitational binding energy Eg of the proto-neutron star (of mass 
MpNs), given by 

Eg = _ 3_Q ^ ^q53 

5Rns 

Here Mp^s = I-SMq, Rns = 10 Km and Gn is the Newton’s gravitational constant. 
Out of this 99% of the released energy is carried away by neutrinos, while the rest 1% 
contributes to the kinetic energy of the explosion. To detect this neutrino burst by the 
earth based detector is of primary astrophysical interest of the core-collapse supernova. 
About 10^^ ergs energy was released in the supernova SN1987A explosion in a couple of 
seconds and two collaborations Kamiokande [7j and IMB[S] hrst detected this neutrino flux 
using their earth based detectors. The observed neutrino luminosity in the detector(IMB 
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or Kamiokande) is ~ 3 x 10®^ erg s“^ (including 3 generations of neutrinos and anti¬ 
neutrinos i.e. Vf.-! t'/i, Vt and 14 , z^, 14 ). So Lj, = ^ ~ 3 x 10^^ erg s“^. The mass of a 
typical proto-neutron star Mp^s = ^-^Mq = 3 x 10^^ g. So, the average energy loss per 
unit mass is ~ 1 x 10^® erg g“^s“^. Note that this is the energy carried away by each 

of the above 6 (anti)-neutrino species. Now besides neutrino, if KK graviton, KK radion, 
axion also take away energy, the energy-loss rate due to these new channels enew should be 
less than the above average energy loss rate |5] i.e. 

enew < 10 ^® erg g”^s"^ ( 2 ) 

and this follows from the observed neutrino luminosity per species (total 6 neutrino and anti¬ 
neutrinos, three type each). So, the upper bound on enew (G- Raffelt’s criterion) is a data- 
driven entity and it does not have an implicit dependence on the fact the ensemble of particles 
obey the normal Bose-Einstein(BE) statistics or Fermi-Dirac(FD) Statistics characterized 
by the equilibrium temperature Tsn or by Tsallis statistics where the fluctuation around 
Tsn is taken into account. If any energy-loss mechanism has an emissivity greater than 
IQis erg g“^ s“^, then it will remove sufficient energy from the explosion to invalidate the 
current understanding of core-collapse supernova. 

Using the Raffelt’s criteria of the supernova energy loss rate for any new physics channel, 
we now constrain the scale A of the dark matter effective theory. Now since the core 
temperature of the supernova is fluctuating, we will work here within the formalism of 
g-deformed statistics [9] where this temperature fluctuation is taken into account. 


B. Fluctuating temperature and Tsallis statistics and temperature fluctuation 


The temperature (T) fluctuations in the g-deformed statistics [lO] takes the distribu¬ 
tion of the following form 




1 

rTf) 



exp 



(3) 


where n is the degree of the distribution and /d = ^- The average of the fluctuating inverse 
temperature {3 can be estimated as 


iP) = n{X^) = / f3f{(3)dp = Po 


(4) 
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Taking into account the local temperature fluctuation, integrating over all 13, we And the 
g-generalized relativistic( with particle energy E = + m^c^) Maxwell-Boltzmann dis¬ 

tribution 


V{E) ~ ^(5) 

(1 4- b{q - 1 )E)- 3 -i 

where g = 1 -|- and b = 4 ^- It’s generalization to Fermi-Dirac and Bose-Einstein 
distribution is worked out in m- The average occupation number of any particle within 
this g-deformed statistics ( Tsallis statistics i) formalism, is given by fi{(3,Ei) (z = 1 , 2 
corresponds to particles) where 


^- 

(1 (g - l)6Ei)-3-i ± 1 


( 6 ) 


where the — sign is for bosons and -|- sign is for fermions. Note that the effective 

_ ]_ 

Boltzmann factor x* = (1 -|- (g — l)bEi) ' 3-1 approaches to the ordinary Boltzmann factor 
^-bEi(^_ ^- 0 oEi'^ as g —)■ 1. The g-deformed statistics finds important application in collider 
physics and astrophysics: Beck et a/. [12] use the g-deformed statistics in order to explain 
the measured energy spectrum of primary cosmic rays. With b~^ = kTo = 107 MeV and the 
deformation parameter g = 1.215 and g = 1.222,respectively, they were able to explain quite 
well the flux rate i.e. the upper(upto the knee) portion and the lower (the ankle) portion 
of the cosmic ray spectrum m- Bediaga et abuse the g-deformed statistics to explain the 
differential cross-section for transverse momenta in electron-positron annihilation [T3]. The 
applications of g-deformed statistics for chaotically quantized scalar fields HI, dark energy 
Ha are available in the literature. Recently, Das et a7|T6| found that an ultra-light ra- 
dion in brane-world Randall-Sundrum model produced in the supernova core can also take 
away energy released in SN1987A explosion and thus can explain the SN1987A energy loss 
rate provided g lies within the range 1.18 < g < 1.32. Here we investigate the supernova 
SN1987A cooling due to free streaming of fermionic dark matter. 


III. DARK MATTER PAIR PRODUCTION INSIDE SUPERNOVA 

Electrons are abundant in supernova. The dark matter fermion may be pair produced in 
the s-channel annihilation of electron and positron: e~{pi)e~^{p 2 ) —^ x{P 3 )x{P 4 )■ 
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FIG. 1; Feynman diagram for the process e e"*“ 


XX 


The effective Lagrangian of describing photon( 7 ) and dark matter fermion (x) interaction 
is given by 

^ (7) 

where — d'^A'^, the e.m.held strength tensor. Here /i^ and correspond to 

the magnetic dipole moment and the electric dipole moment of the dark matter fermion y. 

= I[ 7 ^, 7^] is the spin tensor. The 4-momentum vectors of the initial and final state 
particles in the centre-of-mass frame are given by 


Pi = (E,0,0,p^); p2 = (^,0,0,-p^); 

P 3 = {E',psm9cos(j),psm9sm(f),pcos9); 

P4 = {E', —psindcos^, —psin 6 'sin(/), —pcosd). 

The spin-averaged amplitude-square for the process e“(pi)e’''(p2) x(P3)x(P4) is given by 


= 47 rQ;[/i^ {s(l — cos^ 9) + 4m^(l -f cos^ d)} -|- {(s — 4m^)(l — cos^ d)}] ( 8 ) 


The differential cross-section for the process is 


da 


e e 


xx) = 


dVt 

Finally, the total cross-section is given by 


4m^ - 

1 - 


(9) 


cr e e 






[p\{s + 8 m^) d\{s - 4m^)] 


( 10 ) 


Here rriy^ is the dark matter mass, a = ^ and s = (pi +P 2 )^ = (pa +P 4 )^ is the Mandelstam 
variable. 
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Energy Loss rate 


The supernova energy loss rate due to dark matter fermion pair production is given by [5] 

1 


>xx 


' (^e“^e+*^e“e+—^XX ^rel^com) 
PSN 

1 1 


POO /*oo 


PsN 


dEidE 2 


' mv 


EiE2{Ei + -^ 2 )" 

2 D 1 D 2 


^e-e+^XX 


( 11 ) 


where the c.o.m energy Ecom{= Ei + E 2 ) = 2E (where Ei = E 2 = E) and the relative 
velocity V^ei = Psn is the supernova matter density. The cross-section cxe-e+^xx 

and A = (1 + 7 (A - Pi)y + 1 with i = 1,2. Here b = A = ^ 


given in Eq. 
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(we are working in the unit where ks = 1) and r = The electron and positron 

number densities Ug- = f Ei~^ and ne+ = / Introducing the dimensionless 

variables Xi = Ei/T {i = 1,2), we can hnally write the energy loss rate Eq. 11 as, 

aT'^ 


^e-e+^XX 


dxi 


dx 


Xi X2 {xi -f X2) 


7 


127rVs«i:y iiy [(l + ^(Tii-fi,-))" + l] [(1 + + l] 


( 12 ) 


where the function 7 is given by 
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T\xi+X 2 Y 



\^{Xi+X 2 f + 



+ d\ 


|(a:i +X 2 f 



(13) 


Noting the fact that in the g —)■ 1 limit, the g-deformed distribution formula gets converted 
to either the Bose-Einstein or Fermi-Dirac statistical distribution formula (which describes 
the un-deformed scenario) (see the APPENDIX for a proof) i.e. 


Mf3,E,) = 


1^1 


(1 -t- (g — l)bEi) 9-1 ± 1 


.bEi ± X 


g/3oEi X 


(14) 


where with b = = (5o for g —)• 1 and A is the inverse equilibrium tem¬ 

perature To of the supernova core, the energy loss rate in g = 1 case takes the following 
form 


aA ^ ^ Xi{xi+X 2 ) X 2 

[exp (x, - + l] [exp {x, 

where, /re+ = —Pe- 
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IV. NUMERICAL ANALYSIS 


The dark matter produced inside the supernova core via the channel —)■ XX can 

contribute to the supernova energy loss rate, if the emissivity of this channel s(e~e~^ —)■ 
XX) ^ 10^^ crgf g~^s~^ (= 7.288 x 10“^^GeV). 

Since the core temperature(T) of the supernova is fluctuating, we follow the distri¬ 
bution analysis technique |T 0 ] here, where the temperature distribution is characterized by 
it’s mean value T(= T^jy) = 30 MeV (see Section II B for more details about x^ dis¬ 
tribution). Because of this temperature fluctuation, the ensemble of nucleons, electrons, 
dark matter fermions, photons inside the supernova follow a statistics popularly known as 
the g-deformed statistics (or Tsallis statistics [9], see section II for the related discussion), 
which is different than the usual Fermi-Dirac and Bose-Einstein statistics. The parameter 
q characterizing such distribution is called the deformation parameter: for g 7 ^ 1 , it is the 
g-deformed distribution , while for g = 1, it is the regular Fermi-Dirac or Bose-Einstein 
distribution (undeformed distribution). We will investigate here the dependence of the de¬ 
formation parameter g on the scale A of the dark matter effective theory for a dark matter 
fermion of mass lying between 1 — 100 MeV. 

A. Bound on the effective scale A from the e~^ -|- e~ xX process 

Depending on whether the effective coupling of dark matter fermion with photon is char¬ 
acterized by a dipole moment operator of magnetic or electric type, we have the below 
mentioned three cases: 

1. Case I: /ix(~ I/^m) 7^ 0, dx(~ = 0- 

2. Case II: /ix(~ 1/^^) = 0, (i^(~ l/A^) 7 ^ 0. 

3. Case III: ju^(^ 7 ^ 0) '^x(~ 7 ^ 0- Here A^ = A^ = A 

In each case we have two possible scenarios: (i) Scenario I: The deformation parameter 
g > 1 which is known as the g-deformed scenario and (ii) Scenario II: The deformation 
parameter g = 1 which is the undeformed scenario and particles obey the usual Fermi-Dirac 
or Bose-Einstein statistics. 


1. Scenario I: q-deformed statistics (q> 1) 


Using the Raffelt’s criteria e{e —>■ xx) < 10^^ g ^ 7.288 x 10 ^’^GeV) and 


equation (12), we derive a lower bound on the scale A (which we denote as and A^^ in Case 
I and Case II, respectively) of the dark matter effective theory. In Fig. we have plotted 
A^, Arf (Case I, Case II) and A (Case III) as a function of the dark matter fermion mass 
corresponding to g = 1.05 (left hgure) and q = 1.1 (right hgure), respectively. The region 




m;j,(MeV) 


m,(^(MeV) 


FIG. 2: A^_rf fin TeVj (lower two curves) and A. (in TeV) (topmost eurve) are plotted against 
(in MeV) for q = l.t)b(left figure) and q = l.lfright figure), respectively. 


lying on the left of the vertical line corresponding to = 30 MeV and above the horizontal 
curves of each hgure is allowed. For 10 MeV < m^ < 30 MeV [6], the effective scale A 
corresponding to that window (and above the curves) is allowed. Note that A^ ^ remains 
constant (in both hgures) till = 10 MeV and starts decreasing after that. On the left 
hgure (g = 1.05), till = 10 MeV we hnd A^ = A^ = 8.8 x 10® TeV (Case I and Case II) 
and A = 1.24 x 10”^ TeV (Case III). It becomes A^ = 8.7 x 10® TeV and A^ = 8.4 x 10® TeV 
(in Case I and Case II) and A = 1.2 x 10^ TeV (in Case III) at = 30 MeV. As we go 
from left to right hgure, A^^^ gets increased at a given m^. For = 30 MeV, the topmost 
curve gives A = 3.23 x 10^ TeV, whereas the lower two curves give A^ = 2.3 x 10^ TeV and 
Arf = 2.26 X 10^ TeV, respectively. 

In Fig. [^we have shown A as a function of the deformation parameter g corresponding 
to = 5, 10, 20 and 30 MeV for the above three cases (using Rahelt’s criteria and 


analyzing the equation (12)). On the right hgure, the same plots are shown corresponding 
to 1.0 < g < 1.04. The following observations are in order: 
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q 
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FIG. 3: In the left figure, (in TeV) (lower curves) and A (in TeV) (upper curves) are plotted 
as a function of q for = 5, 10, 20 and 30 MeV. On the right, we have plotted the same but for 
1.04 <q<l. 


1. For a given m^, the lower bound on A increases as q increases. For example, for 

= 30 MeV, A changes from 4.8 x 10® TeV to 3.4 x 10^ TeV as q changes from 1.01 
to 1.1 (left hgure). 

2. For a given q, A decreases with the increase in m^. For example, for q = 1.03, we 
see that as increases from 10 MeV to 30 MeV, A decreases from 8 x 10® TeV to 
7.7 X 10® TeV. 


2. Scenario II: Undeformed statistics(q = 1) 

The q-deformed distribution becomes the usual (undeformed) Bose-Einstein or Fermi- 
Dirac type distribution for g = l(in Appendix A, we have given a derivation of this). We 
are now to investigate how the bound on A gets changed as one switches from a deformed 
distribution to an un-deformed distribution. In Fig. we have plotted and A against 
for g = 1 using the Raffelt’s criteria and equation ( [l5| ). The region on the left of the vertical 
line at = 30 MeV and above the horizontal curves is allowed. Up to = 10 MeV, we 
hnd A^ (Case I) and A^ (Case II) are found to be constant at ~ 2.6 x 10® TeV and after 
which they start decreasing. The topmost curve corresponds to A (Case III) for different 
values of m^. For example, corresponding to = 1 — 10 MeV, we hnd A ~ 3.6 X 10® TeV, 
while at = 30 MeV, we hnd A = 3.3 x 10® TeV. We have summarized our result in Table 
1. From Table 1, we see that as increases, the bound on A^ and A decreases. 
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FIG. 4: Af^^(i(in TeV) and A(in TeV) are plotted against (in MeV) for q = 1. 


Table 1 


mx (MeV) 

Case I : A^ (TeV) 

Case II ; A^ (TeV) 

Case III : A (TeV) 

10 

2.6 X 10® 

2.6 X 10® 

3.6 X 10® 

30 

2.4 X 10® 

2.3 X 10® 

3.3 X 10® 


Table 1: The lower bound on the effective scale A^, A^ and A (TeV) are shown for different 
dark matter mass (MeV) in the undeformed scenario. 

The lower bound obtained on is comparable with that obtained by Kadota et al. mi- 

V. DARK PARTICLE FREE STREAMING/TRAPPING 

The constraint on A obtained in earlier section holds to be true if the produced dark 
matter fermion free streams out of the supernova. To hnd the free streaming let us calculate 
their mean free path |T8] 

iLe 

where ne{= 8.7 x is the number density of the colliding electrons in the supernova 

and aex^ex is the cross section for the scattering of dark matter fermion on electron which 
is related via the crossing symmetry to the annihilation cross-section aee^xx- Now, most 
of the dark matter particles produced in the outermost 10% of the star (0.9i?c < r < Rf) 
from electron-positron annihilation (TH]. Then any of the dark matter particles produced 
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in electron-positron annihilation while propagating through the proto-neutron star, can un¬ 
dergo scattering due the presence of neutrons and electrons inside the star. In the case of 
supernova cooling, neutron-dark matter particle scattering will be negligible for free stream¬ 
ing due to neutron mass [IE]. We use the optical depth criteria [IH] 

dr 2 
^"3 

to investigate whether the dark matter fermion produced at a depth free streams out of 
the supernova and takes away the released energy or getting trapped inside the supernova. 
Here we set ro = 0.9i?c in our analysis, where Rc{— 10 km) is the radius of the supernova 
core (proto-neutron star) [18] . From the optical depth criteria, we hnd that minimum length 
of the mean free path for free streaming A/s is A™"’ = 1.5 km and it increases with A. In 





FIG. 5: On the left side, the free streaming length A/s (km) is plotted as a function of A (GeV) 
for different m^. We have taken Rc{= R) = 10 km. On the right side, the free streaming length 
A/s (km) is plotted as a function of (MeV). The trapping and free streaming regions are shown 
separately. 

Fig. I^on the left side), we have shown the free streaming length A/s (in km) as a function 
of A for different m^. Results for three different cases: Case I; A,, = A, W = 0 and 
Case II; A^ = A, = 0 and Case III; A^ = A^ = A are shown. The horizontal line 
correspond to A™” = 1.5 km and it’s intersection with the curves gives the lower bound 
on A^(Arf) = 7.42 x 10^ TeV, in Case I (Case II) and A = 1.05 x 10® TeV in Case III, 
respectively. On the right side (of Figure]^, we have plotted the free streaming length A/s 
against m^. The horizontal lines corresponding to A™” = 1.5 km remain constant until 
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= 30 MeV, after which they start decreasing except the Case-II in which it remains to 
be constant beyond = 30 MeV. The region below the line = 1.5 km corresponds to 
the trapping region, whereas the region above this horizontal line corresponds to the free 
streaming region. In Case II, where the dark matter-electron coupling is due to electric 
dipole moment of the dark matter fermion, the free streaming of dark matter particles is 
allowed for the entire mass range and thus are not trapped inside the supernova. In Fig. 

we have plotted the lower bound on A against the dark matter mass (considering 
Tsn = 30 MeV, psN = 3 x 10^"^ gm/cc, Ue = 8.7 x m“^). 


2 X 10“ 

1 X 10* 

5 X lO’ 

> 2X10’ 

f- 

^ lx 10’ 

5 X 10‘ 

2 X 10“ 

1 X 10‘ 

FIG. 6: A (in TeV) is plotted against (MeV). The lower set of curves follow from the SN1987A 
cooling in undeformed scenario(q=1.0), while the upper set of curves follows from the free streaming 
of dark matter fermions from the supernova core. The middle set of curves follows from the 
SN1987A cooling in q-deformed scenario(q=l.l). 

The upper(most) set of curves corresponds to that obtained using optical depth criteria 
(based on free streaming of dark matter fermion from the supernova core), whereas the 
lower(most) set of curves are obtained after applying the Raffelt’s criteria on the SN1987A 
energy loss rate in undefromed case and the set of curves in the middle are obtained from 
the SN1987A cooling in q-deformed scenario(q=l.l). The region above both the curves 
(free streaming and SN cooling) is allowed for two different scenarios, q=1.0 and q=l.l. 
correspond to = 10 MeV and 30 MeV. For a dark matter of mass rriy^ = 30 MeV, the 
SN1987A energy loss rate gives a lower bound A(= A^) = 2.4 x 10® TeV (lower curve) in 
undeformed scenario(q=1.0) and A(= A^) = 2.3 x 10^ TeV (lower curve) in q-deformed 
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scenario(q=l.l), whereas optical depth criterion (i.e. free streaming) hxes the lower bound 
on A(= A^) at 7.42 x 10’^ TeV (lower curve). 


VI. DEPENDENCY OF THE EFFECTIVE SCALE A ON THE SUPERNOVA 
PROPERTIES 

So far in our discussions we have not considered the variation of supernova properties like 
temperature(T sat), matter density(p 5 jv), number density of electron(ne) which are relevant 
in the ey —)■ ey scattering process and on the variation of the lower bound on the effective 
scale A. The analysis was based on the following assumptions 

Tsn = 30 MeV, psN = 3 x 10^"^ gm/cc, Ue = 8.7 x 10^^ m“^ 

where, we considered Tsn is the average temperature of the supernova, psN is the average 
matter density of the supernova and Ue is the number density of the colliding electrons in 
the supernova (which are taking part in the ey —)■ ey scattering process). 

In the analysis below, we propose the following 

(i) We can look for the fermionic DM emission pattern which is consistent with the 
Kamiokande [7] and IMB [S] data for neutrino emission due to supernova SN1987A 
explosion (which is more relevant as per our current understanding of supernova ex¬ 
plosion mechanism). 

(ii) In another possibility, we can think about the annihilation process —>■ yy which is 

happening only at the outer 10% of the proto-neutron star at some higher temperature 
(say 50 MeV) and some lower density (say psN = 10^^ gm/cc) [T91 - I23] . The justihcation 
of this second possible assumption is due to infall of matter particles in the accretion 
phase of the supernova, temperature is increasing at the outermost (10%) region of the 
supernova core. Also due to the electron capture process by the protons, the central 
region of the proto-neutron star will be neutron-rich region and will not have much 
free electron-positron pair which can take part in the annihilation process. 
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A. Case I 


We can consider the supernova cooling phenomena due to the free streaming of dark 
matter fermions in both the phases, accretion phase and Kelvin-Helmholtz cooling phase. 
Due to infall of the matter objects in the accretion phase the temperature at the outermost 
region will increase than the average. The number density(ne) of the colliding electron 
(which takes part in ex —t ex scattering process) although appears to increase due to 
infall but actually not so. The infall boosts the production process of dark fermions(due to 
electron-positron annihilation: —)■ XX ) the electron capture process by protons to 

create a neutron-rich core- together (annihilation and capture processes) result the number 
density(rze) to fall to a lower value. 

In Fig. we have plotted the electron number density Ue (in m~^) against the mean 
free path A/^ (in km) for q = 1.0(undeformed scenario) and q = 1.1 (deformed scenario), 
respectively. 
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FIG. 7: The electron number density Uf, (in m~^) is plotted against the electron’s mean free path 
Xfs (in km). We have taken the dark matter fermion of mass m^. = 30 MeV and the value of the 
deformation parameter q = 1.0 (lower curve) and q = 1.1 (upper curve), respectively. 

We have considered the average supernova temperature T = 30 MeV and the number 
density Ue = 8.7 x 10^^ m“^ [THj. Considering the fact that the dark matter fermions, 
produced at a depth rg = 0.9Rc, free stream out of the supernova and take away the 
released energy, we hnd (from the optical depth criterion |T8]) that the minimum length of 
the mean free path for free streaming (A/^) is A™"" = 1.5 km. 
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From figure [^, we see that if the value of the mean free path ~ 1.5 km, we need to 
choose Ue ~ 10^^ m“^ in the undeformed scenario(q=1.0) and Ue ~ 10^^ m“^ in the q- 
deformed scenario(q=l.l), respectively. Suppose at the outermost region the number density 
drops to ~ 10^^ m“^ (from the average value 8.7 x 10^^ m“^ ) and temperature rises to 
T = 50 MeV (above the average value T = 30 MeV ) which is consistent with existing 
supernova simulation (matter density also decreases from its average value 3 x 10^^ gm/cc 
to 10^^ gm/cc) [20] (we compared with the s23WH07 model as in that model mass of the 
progenitor was considered to be 23Mq and progenitor of SN1987A had mass ~ 2 OM 0 ). In 
Fig. we have plotted the lower bound on A against the dark matter mass (obtained 
from the supernova energy loss rate and from the free streaming of the produced dark 
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> 

lx 10 ’ 
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FIG. 8: The lower bound on A (in TeV) (obtained from energy loss rate and the free streaming 
eriteria) is plotted against (in MeV). The lower set of curves follows from the SN1987A cooling 
for q=1.0 and Tsn = 30 MeV, while the curves in the above follow from the free streaming of 
dark matter fermions from the outermost region of the supernova core. The SNe cooling curves 
in q-deformed scenario(q=1.1) (the curves in the middle) are almost overlapping with the curves 
followed from free streaming. 

fermions from the outermost region of the supernova core). The free streaming process is 
insensitive to the temperature Tsn, we set psN = 3 x 10^^ gm/cc and Ue = 1 x 10^^ m“^ in 
the free streaming case. 

The upper set of curves are obtained from the optical depth criteria (based on free 
streaming of dark matter fermion from the outermost part of the supernova core), whereas 
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the lower set of curves are obtained by applying the Raffelt’s criteria on the SN1987A 
energy loss rate. The free streaming curves have shifted downwards from the one obtained 
in £gure|^ and the SN cooling curves in q-deformed scenario(q=l.l) is almost merging 
with the set of curves for free streaming. The two vertical lines correspond to = 10 
MeV and 30 MeV, respectively. For a dark matter of mass = 30 MeV, the SN1987A 
energy loss rate gives a lower bound A(= A^ = A^) = 3.3 x 10® TeV (upper curve) in 
the undeformed scenario(q=1.0) and A(= A^ = A^^) = 3.2 x 10^ TeV (upper curve) in the 
deformed scenario(q=l.l), whereas the optical depth criterion (i.e. free streaming) hxes the 
lower bound on A(= A^ = A^^) at 3.6 x 10^ TeV (upper curve). This suggests almost all 
of the produced dark matter fermions can freely stream out of supernova to contribute to 
the supernova cooling with the fact that they are obeying q-deformed statistics: q=l.l and 
with Ue = 1 X 10^^ m“^), which was not the case for T = 30 MeV in undeformed scenario. 

For the crust temperature T = 50 MeV and the number density of colliding electron 
Ug = lO'^^ m“^ we have plotted the mean free path(A/g) as a function of the effective scale(A) 
and temperature(T) in hgurej^ for different values of dark matter fermion mass m^. 




FIG. 9; In the left figure Xfs(in km) is plotted against A(in TeV) for T = 50 MeV whereas in the 
right figure X(in km) is plotted against T(in MeV) for = A = 3.23 x lO'^ TeV (bound 

obtained from the SN cooling case in q-deformed scenario with q=l.l). 

In £gure|^(left) we see that the mean free path Xfs increases with the effective scale A 
for different dark matter fermion mass. On the right, we see that it hrst increases with 
temperature but eventually becomes constant with temperature T and they are different for 
different (GeV). 

We next study the variation of the effective scale A with temperature for different mass of 
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dark matter fermions for a given A js and different values of mean free path for a given rriy^ and 
are shown in hgure [lo] . From £gure[^ it is clear that the effective scale A varies feebly with 




FIG. 10: A = = Ad (in TeV) is plotted against T (in MeV) where we have chosen = 

10^3 m 3. The left figure corresponds to the variation of A with respect to T for fixed mean free 
path X = 1.5 km and for different whereas the right figure shows the variation of A with respect 
to T for different mean free path X (in km) and for fixed m^. = 30 MeV. 

temperature, it remains almost constant. Variation of the effective scale A with the number 
density of colliding electrons ne (which are taking part in the ex —t ey scattering process) has 
been shown in hgure [TI] . We find that A increases considerably with the increment of the 



FIG. 11: The effective scale A = A^ = Ad (in TeV) obtained for free streaming cases (using 
egn.) %s plotted against rtf. (in. in ) for Ay^ — 1.5 km^ — 30 AfclA We find A. 4 x 10^ 
TlsF for Ue = 10^3 ixi“3 and A ~ 1 x 10® TeV for Ue = 8.7 x 10^3 ni“3. 

number density of colliding electrons tif. (which are relevant in the discussion of the ey —)■ ey 
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scattering process). For rzg = 10^^ m we find A ~ 4 x 10^ TeV and for Ue = 8.7 x 10^^ m 
A ~ 1 X 10® TeV. 


B. Case II 


In section V, we considered the supernova core temperature Tsn = 30 MeV and the 
core density as psN = 3 x 10^"^ gm/cc. From the equations (12, 13) it is clear that, the 


lower bound on the effective scale A depends on the properties (temperature and density) 
of the supernova. As free streaming is happening from the outermost 10% of the supernova 
core, so it is justihed to consider only crust (i.e., crust temperature and density) in the 
discussion of the supernova cooling phenomena as well. For core-collapse supernova, the 
crust temperature will be higher than the core whereas the density will fall from core to 
crust [201 IMl I2S]- In this section, we calculate the lower bound on A in un-deformed (g = 1) 


scenario using Raffelt’s criteria and equation (15) considering temperature T = 50 MeV 
and density p = 10^^ gm/cc at the outermost sector (i.e. at r = 0.9 i?c, the crust) of the 
supernova core. 
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FIG. 12; The lower bound A = A^j, = (TeV) (in TeV) is plotted against the dark matter fermion 
mass (MeV) for two cases: (i) free streaming of dark matter fermions from the crust and (ii) 
from the Raffelt’s criteria of SN1987A energy loss rate. 


In Fig. 


12 


we plot the lower bound on A(= A^ = A^) (TeV) as a function of 


The almost parallel set of curves (also parallel to x- axis) obtained using optical depth 
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criteria (based on the free streaming of dark matter fermion from the ontermost sector 
of the supernova core), whereas the other set of curves are obtained after applying the 
Raffelt’s criteria on the SN1987A energy loss rate (equationj^). The region above each 
curve is allowed. The two vertical lines correspond to = 10 MeV and 30 MeV. For 
a dark matter of mass = 30 MeV, the SN1987A energy loss rate gives a lower bound 
A(= A^ = Ad) = 1.42 X 10® TeV (upper curve), whereas optical depth criterion (i.e. free 
streaming) hxes the lower bound on A(= A^ = A^) at 1.05 x 10® TeV (upper curve). Note 
that the dark matter fermions produced in the outermost sector(i.e. at 0.9i?c) can freely 
stream out to contribute in the supernova cooling phenomena which is not the case if they 
are produced at some inner region (i.e. at distance r < 0.9i?c) of the supernova core. Even 
they are be copiously produced, they are not allowed to free stream (restricted by optical 
depth criterion) and hence can’t contribute to the SN cooling as is seen from Figure]^. Some 
of them will transfer part of its energy back to the medium via scattering and as a result, 
the cooling process occurs slowly. 

Below in Table 2 and Table 3, we summarize our result and make a comparative study 
between two cases as discussed in section V and section VI: 


Table 2 


SN Properties 

(MeV) 

A = A^ = Ad (TeV) 

Free Streaming 

SN Cooling 

Tsd\f — 30 MeV, 

p = 3 X 10^^ gm/cc 

10 

1.05 X 10® 

3.66 X 10® 

30 

1.05 X 10® 

3.34 X 10® 

Ts]\f = 50 MeV, 

p = 10^“^ gm/cc 

10 

1.05 X 10® 

1.45 X 10® 

30 

1.05 X 10® 

1.42 X 10® 


Table 2: The lower bound on the effective scale A = A^ = A^ (TeV) (obtained from Rajfelt’s 
criteria and optical depth criterion) are shown for different dark matter mass (MeV) 
in the undeformed scenario for different values of SN temperature T(= Tsn) and matter 
density psN- 

From Table 2, we see that as increases from 10 MeV to 30 MeV, the bound on A 
decreases (follows from SN energy loss rate) whereas for free streaming case (optical depth 
criterion) A remains constant. We also see from Table 2 that the bound on A change as a 
function of the supernova properties (i.e. it’s temperature and density). In Table 3, we have 
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shown the lower bound A = (TeV) (obtained from Raffelt’s criteria and optical 

depth criterion) for different dark matter mass (MeV) in the undeformed scenario and q- 
deformed scenario for different values of supernova crust temperature T and different values 
of the electron density Ue- 


Table 3 


SN Properties 

(MeV) 

A = A^ = Aa (TeV) 

Free Streaming 

SN Cooling(q=1.0) 

SN Cooling(q=l.l) 

T = 30 MeV, 

Ue = 8.7 X 10"^^ m“^ 

10 

1.05 X 10® 

3.66 X 10® 

3.23 X 10^ 

30 

1.05 X 10® 

3.34 X 10® 

3.23 X 10^ 

T = 50 MeV, 

TT-e = 1 X 10"^^ m“^ 

10 

3.6 X 10^ 

3.66 X 10® 

3.23 X 10^ 

30 

3.6 X 10^ 

3.34 X 10® 

3.23 X 10^ 


Table 3: The lower bound on the ejfeetive scale A = A^ = A^ (TeV) (obtained from Raffelt’s 
criteria and optical depth criterion) are shown for different dark matter mass (MeV) in 
the undeformed scenario and q-deformed scenario for different values of SN crust temperature 
T. We have considered the variation of the temperature and the number density of colliding 
electrons Ue only at the crust level. 

If we consider different phases of a typical core-collapse supernova like the accretion and 
Kelvin-Helmholtz cooling phase, then the bounds will change. We can compare the bounds 
with the knowledge of temperature and density of supernova in those phases. 


VII. CONCLUSION 

The dark matter fermions, pair produced in electron-positron collision e’''e“ —)■ XX inside 
the supernova core, can take away the energy released in the supernova SN1987A explosion. 
Working within the formalism of g-deformed statistics (as the core supernovae temperature 
is fluctuating with the average value Tsn = 30 MeV) and using the Raffelt’s criterion on 
the emissivity for any new channel e{e~^e~ —)■ xx) ^ 10^^ erg we hnd that as the 

deformation parameter q changes from 1.0 (undeformed scenario) to 1.1 (deformed scenario), 
the lower bound on the scale A of the dark matter effective theory varies from 3.3 x 10® 
TeV to 3.2 X 10^ TeV for a dark matter fermion of mass = 30 MeV. Using the optical 
depth criteria on the free streaming of the dark matter fermion, we hnd the lower bound 
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on A ~ 10® TeV for = 30 MeV. In a scenerio,where the dark matter fermions are 
pair produced in electron-positron annihilation in the outermost sector of the supernova 
core (with radius 0.9i?c < < -Rc where Rc{= 10 km) being the supernova core radius or 

the radius of proto-neutron star), we hnd that the bound on A (~ 3 X 10^ TeV) obtained 
from SN cooling criteria (Raffelt’s criteria) is comparable with the bound obtained from 
free streaming (optical depth criterion) for light fermion dark matter of mass = 10 — 30 
MeV. In a nutshell, all the dark matter fermions produced in the outermost sector (i.e. at 
the crust) can freely stream out to contribute to the supernova cooling phenomena which is 
not the case if they are producing in some inner region than the crust. 
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IX. FEYNMAN RULES 

Process: e“e+ —^ XX'- 
e“ e"*" —)■ 7 vertex: 

1 ^ XX vertex: i + dxO'^'^QuX^) 

X. APPENDIX: 

Prom q-deformed statistics to undeformed scenario 

In general, the distribution function for the q-deformed statistics is [9] 

Di = “I" ~ “I" 1 (1®) 

with b = 4 ^ 1 ^, /So = (we work in the unit ks = 1) and r = ^ 

In terms of the dimensionless quantity x* = ^ 

Di = {l + h{q-l){Txi-+ l (19) 
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Now replacing g — 1 by m, (m —)■ 0 as g —)■ 1) 


(1 + b{q - l){Txi - = (1 + bm{Txi - /i*))- = |/(say) 


( 20 ) 


Now 


lim y 

m-5>0 

lim Inu 

m -?>0 


lim (1 + bm{Txi — /i*))™ 

m—>-0 

lim — In (1 + bmiTxi — Hi)) 
m 


= lim --—T-- 

1 + bm[Txi — Hi) 

= b(Txi - Hi) 


b{Txi - Hi) 


Also for g —)■ 1, we find b{= ■^^) = /So = So we find 


lim In?/ = l3o{Txi - Hi) 

q^l 


lim?/ = exp 

q-^l 

= exp 

Clearly, in the undeformed scenario (i.e. g = 1) 


-{Txi- Hi) 


Hi 

Xi — — 


T 


limD, = 

<7^1 


lim y + 1 

q^l 


= exp 


^i- Y 


+ 1 [Proved] 
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